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Experimental bounds on squarks of the first two generations assume their masses to be eightfold 
degenerate, and consequently constrain them to be heavier than ~ 1.4 TeV when the gluino is lighter 
than 2.5 TeV. The assumption of squark-mass universality is neither a direct consequence of Minimal 
Flavor Violation (MFV), which allows for splittings within squark generations, nor a prediction of 
supersymmetric alignment models, which allow for splittings between generations. We reinterpret a 
recent CMS multijet plus missing energy search allowing for deviations from (7(2) universality, and 
find significantly weakened squark bounds: a 400 GeV second-generation squark singlet is allowed, 
even with exclusive decays to a massless neutralino; and in an MFV scenario, the down-type squark 
singlets can be as light as 600 GeV provided the up-type singlets are pushed up to 1.8 TeV, for a 
1.5 TeV gluino and decoupled doublet squarks. 



Introduction: As a solution to the electroweak gauge 
hierarchy problem in the Standard Model (SM), super- 
symmetry (SUSY) is an immensely compelling paradigm. 
However its most popular incarnation, the Minimal Su- 
persymmetric Standard Model (MSSM), has a huge pa- 
rameter space with over a hundred independent param- 
eters. A rather small subset of this parameter space is 
constrained by naturalness [IHS], and it is possible to 
reinterpret current searches in the context of an effec- 
tive theory containing only the most relevant degrees of 
freedom [HIS]. The requirement of naturalness, however, 
gives no guidance regarding the vast majority of param- 
eters of the MSSM or its extensions: those related to 
the squarks of the first two generations are largely in- 
significant to the naturalness argument due to the small- 
ness of their Yukawa couplings. At face value, ATLAS 
and CMS simplified model searches disfavor first- and 
second-generation squarks below ^ 1.4 TeV, for a gluino 
mass of 2.5 TeV or less. However, these limits assume 
an eightfold degeneracy for the masses of two flavors of 
electroweak doublets, Qi,2, and up- and down- type sin- 
glets, (ur, Cfi) and {dji, sn). We will argue below that 
this assumption is not justified. 

One might wonder how drastic the practical conse- 
quences of relaxing this assumption would be. A naive 
rescaling, assuming that the squark production cross 
section goes like m^^, would suggest that the bound 
on a single squark degree of freedom should be around 
8~^/^ ^ 30% smaller than the current limit, hardly a dra- 
matic change. This estimate fails to take into account two 
important effects: the first is the drop in signal efficiency 
at low squark mass due to the hard cuts necessary to 
minimize SM backgrounds; the second is due to Parton 
Distribution Functions (PDFs), which, for non-decoupled 
gluinos, result in the cross sections for second-generation 
squark production being smaller than the correspond- 



ing ones for the first generation (previously noted in the 
context of Dirac gluinos by [H [7]). Both these effects 
work in concert to further weaken the bound on a second- 
generation squark with respect to the naive expectation, 
and make, in themselves, a strong case for reanalyzing 
the data in the context of non-degenerate light squarks. 

Adding the theoretical perspective only makes the case 
more compelling. Squark-mass universality is motivated 
solely by high scale SUSY breaking models such as Mini- 
mal Supergravity, and is not required to solve the SUSY 
flavor problem. Indeed it is not a direct consequence of 
Minimal Flavor Violation (MFV) (see e.g. [3 [9]), which 
allows for splittings between squarks belonging to differ- 
ent representations of the SM gauge group, lifting the 
eightfold degeneracy to a 4 -I- 2 -I- 2 pattern. Further- 
more, SUSY alignment models, which address the SM 
flavor puzzle as well as the SUSY flavor problem, natu- 
rally predict an anarchic sfermion spectrum, with 0(1) 
splittings between squarks. This is achieved by assign- 
ing particular U{1) charges to the different squark gen- 
erations under a new set of flavor ('horizontal') symme- 
tries [To], which forces the squark doublet soft masses to 
'align with' (be diagonal in) the down Yukawa mass ba- 
sis, trivially satisfying the most severe constraints from 
CP violation in the Kaon and B systems and relaxing the 
constraint on the mass splitting (see e.g. [H]). The latter 
was previously thought to be strongly constrained at the 
percent level [121 US], due to a combination of bounds 
from K — K and D — D mixing. However, it was recently 
understood that the bounds from CP-violating processes 
are largely avoided by successful alignment models [14j . 
leading to a much weaker constraint on squark mass non- 
degeneracies, for doublets as well as singlets. 

To summarize, there is strong practical, as well as 
theoretical, motivation for re-examining the experimen- 
tal bounds on the masses of the first two generations of 
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FIG. 1: (Left panel) Variation of signal efficiency x acceptance (e) with squark mass for a single channel in the CMS multijet 
+ MET search with 4.98 fb~^ of data, for a squark simplified model with a 50 GeV neutralino. (Right panel) Cross section 
o (dotted lines), and cr x e for the chosen channel (solid lines) for 8 squark degrees of freedom (in red), naively rescaled for a 
single degree of freedom (in blue). The fiducial cross section limit is (T95 (dotted green). The steeply falling efficiencies at low 
squark masses result in a significant reduction of the limit for the rescaled cross section. 



squarks when the assumption of a full eightfold degener- 
acy is relaxed. In the remainder of this paper we describe 
in detail our methods for reinterpreting current LHC 
searches within this more general framework, and present 
estimated limits for a number of distinct, phenomenolog- 
ically interesting scenarios with non-degenerate squarks. 
We end with a comparison of the effectiveness of the 
different light squark search strategies at ATLAS and 
CMS, and comment on the possibility of optimizing such 
searches by readjusting analysis cuts. 

Procedure: We determine limits on simplified spectra 
consisting of a gluino and one or more squarks, decaying 
to a massless neutralino plus jets, with other superpart- 
ners decoupled. We focus on the recent CMS 7 TeV mul- 
tijet plus missing energy (MET) search with 4.98 fb^^ 
of integrated luminosity [15j : we compare this with the 
performance of other CMS and ATLAS jets plus MET 
searches below. For the squark-neutralino and gluino- 
neutralino simplified models, CMS makes available effi- 
ciency maps in the (niqjm^g), (mg,m-^g) planes: we use 
this information directly whenever possible. Our anal- 
ysis also requires efficiency maps for a squark-gluino- 
neutralino simplified model, with on-shell intermediates 
in decays where relevant, as well as for mixed produc- 
tion of two squarks of different masses. Lacking the 
pertinent efficiency and acceptance information in these 
cases, we simulate them using Pythia 6.4.24 [l^ with 
the CTEQ6L1 PDF set [H], subsequently passing the 
events through two different analysis pipelines, ATOM 
("Automatic Test of Models") [IB], and the PCS [H] 
'theorist-level' detector simulation, as an internal cross- 
check. The former (soon-to-be publicly available) is a 
RIVET ^-based tool that estimates current LHC limits 
on a given model, flagging problematic regions e.g. where 
the signal leaks into control regions, or where signal ef- 
ficiencies are too sensitive to cut positions. We validate 



both pipelines by reproducing the published limits of the 
search within ^ 50 GeV, corresponding to an accuracy 
in the estimated experimental acceptances of better than 
20%. We correct the various production cross sections 
using NLO-t-NLL K-factors from NLLfast [1IH21] in all 
figures except Fig. [2] where we use Prospino 2.1 |21) . 
For mixed production, which is not considered in NLL- 
fast/Prospino, we cannot use such a procedure reliably. 
In this case we present the result using the LO cross sec- 
tion, and show the effect of including K-factors of 1.5 and 
2.0. We have checked using MadGOLEM [,25j that, for a 
selection of points close to the current limit, the chosen 
K-factors are in fact conservative. Finally, at each point 
in signal parameter space, we combine the fourteen non- 
overlapping signal regions using an approximated likeli- 
hood C = niiipoiss(nf^|s, + hi) ■ gauss(6,|6f P, 5b,). All 
limits are 95% confidence-level exclusions derived from 
a profile likelihood ratio using the CLg technique |26) . 
computed with a private toy MC code, verified by com- 
paring it to asymptotic distributions ^27j calculated by 
RooStats [2H]- 

Note that the NLO squark production cross section, 
as computed by Prospino 2.1 in the limit of large gluino 
mass, decouples much more slowly than one would expect 
(cr oc w^^), leading to an overestimate of the K-factors 
in this corner of parameter space. When using Prospino, 
we are careful to stay in the safe region. Note also that 
the naive use of Pythia, with all the squarks and gluino 
production channels enabled and non-degenerate squark 
spectra can result in under-sampling of certain regions of 
phase space. As a workaround, we performed multiple 
runs corresponding to different sub-processes. 

Results: As mentioned above, naively rescaling the 
squark cross section limit by the number of squark de- 
grees of freedom in the first two generations significantly 
underestimates the change in the limit for a single squark 
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FIG. 2; The limits in the squark-gluino plane for production and decay of a single squark degree of freedom of different flavors. 
The effect of the PDF enhancement on the production of a first-generation squark, via a t-channel gluino, is clearly manifest 
in the form of a more stringent limit for fixed gluino mass; this difference approaches zero in the gluino decoupling limit. 



for two non-trivial reasons. The first is the sharp drop in 
experimental efficiencies at small squark masses. This is 
due to hard cuts (mainly on variables correlated with 
the visible and invisible energy of the event, such as 
Meff , Ht , ^ , I^t) placed on the data in order to sup- 
press the large SM backgrounds. We illustrate this point 
in Fig. [l] The left panel shows the variation of the CMS- 
provided efficiency x acceptance (e) with squark mass for 
a chosen channel in the squark simplified model analysis, 
with a neutralino mass of 50 GeV. The steep decrease in 
signal efficiency for squark masses below about 500 GeV 
seen in this plot is typical of all pertinent cut-and-count 
searches, making these particularly inefficient for light 
squarks. The right panel shows the cross section a (dot- 
ted lines), and axe (solid lines) for the chosen channel 
for 8 squark degrees of freedom (in red), naively rescaled 
for a single degree of freedom (in blue), with decoupled 
gluinos. The fiducial cross section limit in this region, 
CT95, is indicated by the dotted green line. The rising 
cross section with decreasing squark mass is compensated 
by the steeply falling efficiency, significantly reducing the 
exclusion limit for the rescaled estimate. Note that this 
is not the true bound on a single squark, but simply il- 
lustrative of the repercussions of the falling efficiencies. 

The second important reason is due to PDFs: the 
large valence quark density in the proton leads to a 
squark pair-production cross section that is dominated 
by first-generation squark production with gluinos in the 
t-channel, if the latter are accessible. The current limit 
for non-decoupled gluinos is therefore driven predomi- 
nantly by first-generation up-type squarks, with limits 
on the second-generation being correspondingly weaker. 
We neglect squark mixing for simplicity (see conclusion 



for comment). The PDF effect is shown in Fig. [2] where 
we plot contours of cr/criim, the ratio of the total cross sec- 
tion to the excluded cross section, in the squark-gluino 
plane, for a 50 GeV neutralino, and a single squark de- 
gree of freedom of varying flavor. The estimated limit is 
much weaker for a second-generation squark, where the 
cross section has no PDF enhancement. This difference 
is expected to decrease with increasing gluino mass, as 
the contribution from the i-channel gluino decreases, but 
due to the slow decoupling of the gluino, the asymptotic 
behavior is reached for gluino masses larger than those 
shown here. 

The striking contrast between the first- and second- 
generation squark bounds is due in part to the fact that 
we are taking an artificial limit in which all other squark 
degrees of freedom are formally decoupled. For a more 
generic, anarchic squark spectrum, one would expect con- 
tributions to the total production cross section from all 
squark flavors. In Fig. [3j we show the estimated limits in 
three more 'realistic' and phenomenologically interesting 
scenarios with multiple non-degenerate squark degrees of 
freedom, and a gluino just above the current limit, at 
1.5 TeV. The left-hand panel covers the least constrained 
scenario of a single light second-generation squark split 
from all the other squark flavors, all of which are ac- 
cessible. The middle panel gives the estimated limits in 
an alignment- type scenario with first-generation squarks 
split from the second generation. It illustrates an inter- 
esting interplay between flavor and collider physics, since 
the splitting between the electroweak doublets cannot be 
arbitrarily large due to the combination of constraints 
from K — K and D ~ D mixing, assuming down align- 
ment [13]. Note that the flavor constraints shown in the 
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FIG. 3: Squark mass limits in three phenomenologically interesting scenarios with non-degenerate first- and second-generation 
squarks. The left panel contains the least constrained scenario, with a single second-generation squark flavor split from all others; 
the middle panel corresponds to an alignment-type scenario with first-generation squarks split from the second-generation. The 
shaded blue region is excluded by flavor and CP violation constraints which apply to electroweak doublet squarks only, while 
the singlet spectrum remains completely unconstrained; the right panel corresponds to an MFV-type scenario with split up-type 
and down-type singlets, and doublets formally decoupled. The red dashed (dotted) lines represent the exclusion contour if the 
LO mixed up-down squark production cross section is multiplied by a K-factor of 1.5 (2.0). 



plot include the full dependence on the squark masses, 
crucial when the splitting is large . Although the sin- 
glet squarks are kept degenerate with the corresponding 
doublets for simplicity, their splittings are unconstrained 
by flavor, and they could also be decoupled, resulting 
in weaker LHC bounds (corresponding to the contour 
f/(7iiin ~ 2), with unchanged flavor bounds. The right- 
hand panel contains the limits in an MFV-type scenario, 
with split up-type and down-type singlets, and doublets 
formally decoupled. The red dashed (dotted) lines rep- 
resent the exclusion contour if the LO mixed up-down 
squark production cross section is multiplied by a K- 
factor of 1.5 (2.0). 

The surprisingly weak limits, in particular for squarks 
of the second generation, demonstrate how ineffective 
current searches are for light squarks. Re-optimizing 
the ATLAS 2-6 jets plus MET search using only the 
nies cut is not effective: while the background grows 
like m^jf, the signal grows much more slowly, ensuring 
that decreasing the nice cut makes things worse. It is 
possible that the limits would improve on performing ei- 
ther a full re-optimization including all cut variables, or 
a shape analysis; such a study, however, is beyond the 
scope of this paper. Instead, in Fig. |4j we compare the 
limits for squark cross sections from various 7 TeV AT- 
LAS and CMS jets plus MET searches (which have limits 
for degenerate squarks that are competetive with those 
of recent 8 TeV searches [331 IS])- We find indeed that 
the most stringent bounds come from the more complex 
shape-based analyses, such as the CMS razor search. 

Conclusion: We have argued that a combination of 
reduced efficiencies and suppression due to PDFs leads 
to constraints on non-degenerate squark masses (for the 
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FIG. 4: Comparison between upper limits on squark pair- 
production cross sections with a decoupled gluino and mass- 
less neutralino, from 7 TeV 5fb"^ ATLAS and CMS jets plus 
MET searches [TH [HOHH^ . We use the official experimental 
limits, except for the ATLAS search where we use our esti- 
mate of the limit, simulating the search with ATOM (solid) 
and PGS (dotted). 

first two generations) that are significantly weaker than 
those assuming eightfold degeneracy. For instance, an 
0(400 GeV) squark belonging to the second generation 
can be buried in the LHC jets plus MET data. In the 
above analysis we have neglected for simplicity the effects 
of squark mixing, which could be sizable in alignment 
models. In addition, our reinterpreted limits, while as- 
suming the bino is the lightest SUSY particle (LSP), are 
still applicable for singlino or gravitino LSPs, or when ad- 
ditional electroweak {e.g. higgsinos) and leptonic states 
are present, but do not drastically alter the light squark 
branching ratios. In spite of the dramatic increase of 
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the squark pair-production cross section with decreasing 
squark mass, the corresponding SM backgrounds tend to 
grow even faster, and dedicated searches will be required 
in order to extract light squarks from the data. One pos- 
sible addition to these is charm tagging, which would im- 
prove the sensitivity for light charm squarks. This is al- 
ready implemented in Tevatron i — >■ cx° searches [55115^ . 
which set limits for a single charm squark which are 
weaker than those obtained here. 

The impact of observing light squarks would go beyond 
'mere' discovery of SUSY (a dramatic event by itself), 
since the knowledge of a split two-generation spectrum 
encodes information about SUSY breaking. Moreover, 
discovering a splitting between generations (breaking of 
the U{2) flavor symmetry) would have sensational impli- 
cations, giving us insight on microscopic flavor dynamics 
far beyond the direct reach of the LHC. 
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